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1. Background and Summary 


Aircraft control systems are usually very reliable because of redundancy and multiple 
control surfaces. However, there are rare occasions when potentially disastrous flight control 
system failures do occur. At such times, the use of appropriate modulation of engine thrust 
to stabilize the aircraft may be the only chance of survival for the people aboard. In several 
cases where complete loss of control systems has occurred in multi-engine aircraft, pilots used 
propulsion system to regain limited control of the aircraft with various degrees of success 
(see [1]). 

In order to evaluate the feasibility of using only engine thrust modulation for emergency 
backup flight control, the NASA Dry den Flight Research Center has been conducting a series 
of analytical studies and flight tests on several different types of aircraft in a propulsion 
controlled aircraft (PC A) program. Simulation studies have included B-720, B-727, MD-11, 
C-402, C-17, F-18, and F-15, and flight tests have included B-747, B-777, MD-11, T-39, Lear 
24, F-18, F-15, T-38, and PA-30 (see [1]— [7]). One objective was to determine the degree 
of control available with manual manipulation (open-loop) of the engine throttles. Flight 
tests and simulations soon showed that a closed loop controller could improve the chances 
of making a safe runway landing (see [2]). The major work to date has concentrated on 
three aircraft (F-15, F-18, and the MD-11). Successful landings using PC A controllers were 
performed on the F-15 and MD-11 without the use of control surfaces. 

During the course of the research, some unique challenges have been identified (see [6]). 
Compared to the conventional flight control surfaces, the engines are slow and have limited 
control effectiveness. Hence the ability of the system to promptly respond to aerodynamic 
changes is limited. Consequently, many nonlinear effects, which are easily accommodated by 
a conventional flight control system, become significant issues in the design of an effective 
controller when the engines are used as the only control effectors. A number of nonlinear 
behaviors observed during PC A flight tests, which are not meant to be exhaustive, are 
reviewed in the next section. 
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The flight controllers for the PCA thus far have been designed based on linearized models 
of the aircraft dynamics and the engines. Nonlinear simulation and flight testing are then 
carried out to tune the controller gains for enhanced performance. We realize, however, that 
an alternative is to recognize the nonlinearities in a PCA, and try to incorporate them as 
much as possible in the controller design. Such a controller is likely to be more responsive 
to the inherent nonlinearities in the system, and thus produces more effective control. This 
investigation represents a preliminary attempt in this direction, and part of the results have 
been reported in the pending article Ref. [8]. In particular, a recently developed nonlinear 
predictive control method is introduced as a potential design tool for the PCA controller. 
This approach bears some of the similar features that make the Linear Quadratic Regulator 
(LQR) method highly successful for linear systems. Thus it appears to be an appealing tool 
to be employed in the design of a nonlinear PCA controller. A problem of controlling a 
crippled aircraft is used to demonstrate the potential of this method in a highly demanding 
environment such as in PCA applications. In this problem, an airplane that lost the use 
of its rudder is controlled with reasonable satisfation by including the nonlinear dynamic 
coupling effects in the control law design. 

As an evolving step toward fully nonlinear PCA controller design, this research also car- 
ried out, in conjunction with a PCA project conducted in the summer of 1996 at Dryden, 
a PCA longitudinal controller design for the C-17 transport aircraft at a specific flight con- 
dition. At that condition the linearized PCA dynamics exhibits the so-called “Tuck-mode” 
behavior which in this case is an unstable mode. The controller that works for other flight 
conditions could not control the aircraft satisfactorily. A successful PCA controller was .Re- 
signed for the C-17 at this condition by an innovative two-step design process. First an 
optimal PID full-state feedback control law was designed to stabilized the aircraft. Then 
this PID control law was used as the starting point for iterative design of an output feedback 
controller using the Nonlinear Control Design (NCD) toolbox in the MATLAB. 1 Some of 
X MATLAB is a trade mark of the Mathworks Inc.. 
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the details are given in Section 5. 


2. Nonlinearities in Propulsion Controlled Aircraft 

The flight control of the PCA is achieved by using differential thrust for lateral/directional 
control and symmetric thrust for longitudinal control. Refs. [2] and [4] give detailed account 
of this principle. Because of the limited control effects of the engines, however, a number of 
nonlinearities in the PCA which are otherwise insignificant when the aerodynamic control 
surfaces are operative become prominent. Some of the notable nonlinear phenomena in the 
PCA include: 

1. Engine dynamics: Engine deceleration inputs usually reduce the thrust faster than 
engine acceleration inputs increase the thrust. The engines are slower to respond when 
the throttles are moved forward than when the throttles are moved back. Typically, 
a high bypass turbofan engine exhibits very slow thrust increase initially. Once the 
the thrust reaches about 20% of the maximum thrust, the thrust response improves 
dramatically. The thrust decay shows similar nonlinear effects. 

2. Engine saturation: Because of the limited control authority, the control commands can 
easily cause the engines to operate at their limits. For instance, a large roll command 
may cause one engine to saturate while the other engine still has more control authority 
available. This will lead to unsatisfactory lateral/directional control. Also, in low-speed 
landing approaches, the commanded engine thrust is close to idle. Disturbances caused 
by gusty conditions could result in engine saturation at the idle. Being a low bandwidth 
system, the engines may also experience rate saturation, and this is compounded by 
the nonlinear aspect of the thrust response. 

3. Ram drag: Ram drag is due to the mass flow rate term which is a function of fan speed, 
nozzle area, exhaust temperature and other flight conditions. It affects the dynamics 
of the PCA and is inherently nonlinear. 
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4. Engine inlet location: This term is aircraft-dependent and has caused the F-15 to 
respond in a minimum-phase fashion; but the F-18 behaves with a nonminimum-phase 
like response depending on the aircraft flight condition. This could in fact be attributed 
to the nonlinear interactions between the airframe and propulsion system. 

5. Dynamic cross-coupling: It has been observed that at high bank angles, the longitu- 
dinal control of the PC A deteriorates dramatically (see [6]). This is caused by the 
cross-coupling effects between the lateral and longitudinal dynamics of the PCA. En- 
gine saturation can also contribute to the cross-coupling which in turn degrades both 
lateral and longitudinal control effectiveness. This degradation of performance results 
from the fact that the longitudinal control lateral control systems are usually designed 
separately based on decoupled linearized dynamics. 

The flight test data of the F-18 System Research Vehicle (SRA) at the NASA Dryden 
Flight Research Center illustrates some of the above aspects. This F-18 SRA, as shown in 
Fig. 1, is a twin-engine, two-seat airplane. The aircraft is powered by two General Electric 
F404-GE-400 afterburning turbofan engines rated at approximately 16,000 lb static thrust at 
sea level. The F-18 features a midwing configuration with a wing-root leading-edge extension 
(LEX) that extends from the forward portion of the fuselage and blends into the wing. 

The F-18 flight test studies were conducted with the normal control surfaces locked in 
a trimmed position and the throttles were given a step input. A series of open loop, pilot 
evaluation maneuvers were flown at a condition of 10,000 ft, V=195 knots, gear up, flaps at 
30 deg., in the direct electric link (DEL) mode with gain override. This mode allowed the 
aircraft to be flown in a ’’true” open loop mode without any feedback parameters interacting 
with engine dynamic responses. 

Figure 2 shows the response to a doublet input in the throttles, denoted by power lever 
angles (PL As). The two PLAs are virtually the same, therefore produce changes in symmetric 
thrust which mostly influences the longitudinal motion. The pitch rate response in Fig. 2 (b) 
decreased first as the PLAs were increased, and then rose up in response to the PLA increase. 
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This would indicate a nonminimum-phase system if this occurred in a linear system. But as 
the example in Section 4 will demonstrate, this could be actually caused by the nonlinear 
PCA dynamics. In Fig. 2 (c) the roll rate showed unsymmetric response to the PLA input: 
when the PLAs were increased, the variation in the roll rate was slow and gentle; when the 
PL As were decreased, however, the roll rate showed a sharp and large transient response. 
This is another indication of nonlinear dynamics of the PCA. 

Figure 3 shows response of the F-18 to a different PLA doublet input. Since the PLAs in 
this case caused differential thrust for the two engines, the lateral /directional response was 
more visible (Fig. 3 (b)), in the roll rate in particular. The engine thrust of the PCA can 
be determined based on several quantities, including fan speed, exhaust gas temperature, 
airflow, discharge pressure and other parameters. The discharge pressures for the two engines 
in this case are shown in Fig. 3 (c). It is clear that the engine response were slower when 
the throttles were increased than when the throttles were decreased. 

3. A Predictive Control Approach 

So far the designs of the controllers for the PCA have been based on linearized aircraft 
and engine models, and linear control methods have been used. In view of the above noted 
significant nonlinear behaviors, it seems logical to expect that an appropriately designed non- 
linear controller may take these nonlinearities into account in the control law, and thus offer 
improved performance for the PCA. Compared with the success of the linear control theory, 
the nonlinear control theory is still in an incomplete and less satisfactory state, although 
remarkable progress has been made in recent years, particularly in geometric control the- 
ory (see [9]) In the aerospace control community, the nonlinear method known as dynamic 
inversion (Ref. [10] and references therein), which is based on the input-output feedback 
linearization technique in geometric control theory, has been the focus of active research. 
But it is felt that dynamic inversion may not be an ideal method for PCA controller design 
because of two major concerns: (1) The control authority of the engines is rather limited. 
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This means that control and control-rate saturations can easily happen. And when they do, 
the control law based on dynamic inversion becomes invalid. (2) The engines, as the control 
effectors, have very complicated nonlinear dynamics, and considerable modeling uncertain- 
ties are bound to be present. Since the principle of dynamic inversion is canceling system 
nonlinearities by state-feedback, robustness could be an issue in the presence of such large 
uncertainties. 

An alternative is a nonlinear predictive control method developed more recently (see 
[11, 12, 13, 14, 15]). The potential advantages of this approach in PCA application lie in 
that no stringent conditions on the system are required, and in the event of saturation, the 
control law remains applicable and the control efforts are distributed in an optimal fashion. 
In the following, the basic concept of this approach and some results are briefly reviewed. 

Consider a multivariable nonlinear control system modeled by 

m 

X = f(x) + G(x)u = f(x) 4 - Yldi( X ) U i (!) 

2=1 

y = c(x) ( 2 ) 

where the state is x , defined in a compact set X C R n , and y R l is the output. A piecewise 

continuous input u £ R m is admissible if for all / > 0 it satisfies the condition 

u(t) ELl = { u{t ) £ R m | Li < Ui(t) < Ui , i = l,...,m} (3) 

where the bounds Li and Ui are given, and can be state- and time-dependent if need be. The 

vector functions f(x ), gi(x), and c(x) are sufficiently differentiable on X. Suppose that a 
desired output is specified by a sufficiently smooth /-dimensional vector function y*(t ), t > 0. 
The requirement on y*(t) is that a corresponding state trajectory x*(t) and control function 
u*(t) exist with #*(/) £ X and u*(t) £ O for all t > 0, such that x*(t) and u*(t) satisfy 
(1), and y*(t ) = c(rr*(/)). While the explicit knowledge of £*(/) and u*(t ) is not needed for 
control law design, their existence ensures that the desired response is feasible. 

The control problem is to design a control law u = r(x,y*,t) G Ll for alH > 0 such that 
for an arbitrary initial conditions of the system (1), denoted by Xq £ Xq C X, y(t ) — >• y*(t) 
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as t — » oo. Note that when y = c(x) = x and y*(t) = x*(t ), the problem becomes a 
state- regulation problem. 

At any t > 0, let 7 ,■ be the relative degree of the ith (1 < % < /) output yi at a?(t) ( 7 ,- is 
then the number of differentiations needed on y^ — Ci(x) until any components of u(t) appear 
explicitly). The influence of the current control u(t) on the system response at t -f h for a 
time increment h > 0 is predicted by expanding each yi(t -f- h) in a 7 ;th order Taylor series 


y(t + h) Pd y(t) + z(x(t), h) + W(x(t), h)u(t) 


(4) 


where each component of z(x, h ) = {z\(x, h ), . . . , z\(x, h)) T is defined by 

h 2 h'n 

Zi — hlPf^Ci) -f — Ty(cj) + -f — LJ (c*), i = l, ..., I 

and the ith row of the / x m matrix W is 

h'yi 

w, = ( — -L gi L' 1 /~ 1 (x i ), — jL gm LJ~ 1 {x i )), i = 1 ,...,/. 
7 *- 7 i- 


(5) 


( 6 ) 


In the above the notation of Lie derivatives has been used to simplify the representation of 
the partial derivatives. Similarly, expand each y*(t + h) in a 7 *th order Taylor series 


y*(t + h) « y*(t) + d(t, h ) 


(7) 


where the ith component of d(t , h ) is 


h 2 

2 ! 


di(t,h) = hy*(t) + '^y*(t) + ... + ^jy* 1 i7i) (t), i = 1 


h* 

7* 


( 8 ) 


The tracking error at the next instant t + h is then predicted as a function of u(t) by 


e(/ + h) = y(t + h) — y*{t + h) « e(t) + h) — d(A h) + W(x(t), u(t ), h). (9) 


The control w(/) is then found by solving the constrained optimization problem 


1 1 
min J = ~e T (t + h)Qe(t + h) + ~u(t) T Ru(t) 

i(t)£Q. 2 2 


( 10 ) 
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where Q and R are positive semidefinite matrices of appropriate dimensions. Define a vector 
saturater s(-) : R m R m by 


Ui, 

Qi > Ui 


& 

II 

s 

Li < qi < U{ , i = 1,2 , ...,m 

(ii) 

Til 

qi < Li 



for any q £ R m . The results regarding the optimal control are summarized in the following 
[11, 12, 13]: 

Replace e(t-\-h ) in (10) by the prediction (9). Suppose that the mxm matrix P(x(t ), h) — 
{ Pij } = W T (x(t), h)QW(x(t), h) + R is nonsingular at x(t). Then, for any h > 0, we have 

1. The unique optimal control u op (t) exists and is the unique solution of the equation 

U = s {pW T Q(d - z - e) - l/3(W T QW + R) - /]«} = p(u) (12) 

where all the quantities are evaluated at t, I is an identity matrix, and 

m m 

/^ = {EEp; j }" I/2 - (13) 

*=1 j = 1 

2. In general the fixed-point iteration sequence { u k } generated by 

u k = p{u k ~ l ), k = 1,2,..., Vn° € R m (14) 

converges to u op (t) globally. 

3. When none of the components of u op {t ) are on the boundary of the control set Ft, the 
closed-form control law from Eq. (12) is 

u = (W t QW + R)~ l [W T Q{d -z- c)]. (15) 

An alternative formulation of the problem can be found in Ref. [14]. Discussions on 
closed-loop stability for several classes of systems under control law (15) and robustness of 
the controller are given in Ref. [12]. 
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The optimal control command at any t > 0 is given by (12) which can be generated by 
the globally rapidly convergent algorithm (14), and this fixed-point algorithm is particularly 
suited for computer implementation. Even in the presence of control saturation, the u(t) 
from the solution of (12) is still optimal in the sense of minimizing the performance index 
(10). When both control and control-rate are bounded, a minor modification in the saturater 
s used in the optimal control law will guarantee strict satisfaction of all the constraints (see 
[15]). When no control saturation is encountered, control law (12) simply reduces to the 
explicit closed form (15). For the PC A control application, this predictive control method 
may prove to be advantageous in this environment with severely limited control-authority 
and significant nonlinearities. 


4. Control of a Crippled Aircraft 


4.1 Aircraft Model 

This example is used to demonstrate the application of the predictive control approach, 
and to illustrate that non-intuitive results may be achieved when the system nonlinearity is 
appropriately taken into account. Consider the following aircraft model which is the same 


as “aircraft A” in Ref. [16]: 

a = —2.9998a + q — (3p (16) 

/? = — 0.05586/? — r + ap (17) 

p = -202.96/? - 39.97p + 2.7943r - 0.70574gr - 601.37$, (18) 

q = -22.692 a - 4.0556,? + 0.71992/?p + 0.95965pr - 61.1294 (T9) 

r = 6.8294/3 — 0.47937r — 0.78674p<? (20) 


where the standard notation is used: a is the angle of attack (deg), j3 sideslip angle (deg), 
p roll rate (deg/sec), q pitch rate (deg/sec), and r the yaw rate (deg/sec). The two controls 
are aileron deflection S a (deg) and elevator deflection S e (deg). 
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This is a nonlinear model with cross-product terms retained, although the control effectors 
are still the conventional control surfaces. Nonetheless, some of the features this model 
presents are quite unconventional because of the nonlinearities, and the notable absence of 
the rudder may be conveniently interpreted as the loss of a control surface (rudder). Hence 
the difficulties in controlling this crippled aircraft resemble to some extent those which would 
be encountered in a PCA. 

First of all, it can be easily shown that the linearized system to Eqs. (16-20) is not 
controllable (the controllability matrix only has rank 3) because the rudder is not available 
for control purpose. So pole-placement or LQR type of techniques for linear controller 
design are not applicable even for the linearized system. In fact, in the linearized system, 
Eqs. (17) and (20) constitute a free stable subsystem in (3 and r which is not influenced by 
any controls and other variables. 2 Secondly, the dynamic inversion approach for nonlinear 
controller design cannot be used for longitudinal control purpose if the usual coordinated 
flight (in which (3 = 0) is desired. This is because with the two controlled longitudinal 
outputs yi — a and t/ 2 = Q, it can be readily verified that the nonlinear system does not 
have a relative degree (see [9]) whenever (3 — 0. 

Since the predictive control method does not require any stringent conditions on the 
system, it appears that this approach is an attractive alternative for flight control design for 
this model of a “crippled” aircraft. We shall consider two control problems in the following. 
The first one demonstrates the effectives of the predictive control method, and the second 
shows that the controller based on the nonlinear dynamics can control the aircraft to achieve 
commanded turning maneuver, which would be impossible to do with only the linearized 
model. 

4.2 Longitudinal Control 

Suppose that the objective is to design a controller to track commands to the longitudinal 
motion while the lateral/directional motion is stabilized at f3 = p = r = 0, assuming full- 
2 This subsystem is nonetheless stable with a light damping of £ = 0.1 and a time constant of 3.9 sec. 
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state feedback. More specifically, suppose that the aircraft is to be controlled to achieve the 
quasi-trim conditions 

a* = cw, q* = q C om , 0* = p* = r* = 0 (21) 

where the commanded constants a com and q com satisfy the relationship 

ow = 0.333355557^. (22) 

This relationship is obtained from setting a — 0. Use the values in Eq. (21) as the reference 
trajectory for a full-state regulation problem ( y = c(x) = x in Eq. (2)). The predictive 
control approach introduced in the preceding section is then employed to design a nonlinear 
controller. The closed-form control law follows directly from Eq. (15). The controller 
parameters are selected to be Q =diag{Qi, ..., Q$}, and R =diag{i?i, ^ 2 } with 

Qi = 1.0, Q 2 = 0.1, Qs = 1.0, Qa = 1.0, ft - 0.1, Ri - 0.1, R 2 = 0.1. (23) 

The value of h = 0.4 is also used. It was found that in this case the system performance 
was insensitive to the choices of R\ and R 2 . Any values of R\ and R 2 in the range of [0, 10] 
would give practically the same response. On the other hand, the choices of Qi seem to have 
more profound influence on the response and stability. A satisfactory set of Qi , however, are 
not difficult to determine by numerical simulations, because these parameters have a similar 
interpretation and effect as the weightings in an LQR problem. 

Figure 4 shows the time histories of the states for the initial conditions 

a(0) = g(0) — 0, P(0) = 5 (deg), p( 0) = 5 (deg/sec), r(0) = —5 (deg/sec) 

and command inputs 

a com = 1.67 (deg), q com = 5 (deg/ sec). 

Figure 5 presents the corresponding control histories. It can be seen that all the lat- 
eral/directional variables subside to zero in about 3 seconds, despite the nonzero initial 
conditions. The longitudinal motion tracks the commanded a com and q com gracefully in 
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about 1.2 second. The local closed-loop stability of the nonlinear system is checked by 
examining the eigenvalues of the linearized system. 

It is interesting to observe the time-history of a(£) in Fig. 4. The initial movement of a(t) 
is opposite to the commanded direction. This would indicate a nonminimum-phase system 
if this phenomenon occurs in the response of a linear system. However, it can be shown that 
the linearized system in this case is actually minimum-phase with all the transmission zeros 
in the left open half of the complex plane. The undershoot of a(t) is caused by a nonlinearity 
in the system, the cross-product term (3p in particular in Eq. (16). Another point is that 
when the predictive controller is constructed based on the linearized dynamics, although 
the linear system is stable, the nonlinear system is not when the controller is applied to the 
nonlinear dynamic equations. Figure 6 shows the state histories of the nonlinear system with 
the linear controller. It is evident that the motion entered a limit cycle. 

The above phenomena clearly indicate that the analysis and controller design based on 
the linearized model in this case are inadequate. 

4.3 Controlled Turn 

Suppose now that the aircraft is to be controlled to achieve a steady-state turn in which 
all the state variables are constant with r = r com . This can be a difficult maneuver for 
this crippled aircraft because the rudder, which is missing, is the primary control surface 
for yaw control. Indeed, if only the linearized dynamics are considered by ignoring the 
nonlinear terms in the dynamic equations, any commanded turn is impossible to achieve for 
any r com ^ 0, since the (3-r subsystem is not influenced by any controls and other states in 
the linearized system. However, when the nonlinearity is retained, this turning maneuver - is 
actually feasible because of the cross coupling. For instance, for r com = 5 (deg/s), such a 
quasi-trim condition is defined by 

a* = 4 (deg), q * = 15.39 (deg/s), (3* = 2.63 (deg), p* = 1.28 (deg/s),r* = 5 (deg/s). (24) 
Use the values in above equation to specify the reference trajectory and the same predictive 
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controller designed for the above pitch control maneuver is employed for flight control. For 
the initial conditions 

a(0) = 13(0) = p{ 0) = q{ 0) = r(0) = 0. 

Figs. 7 and 8 show the histories of the state variables and two controls 8 a and S e . It can 
be seen that the controller effectively drives all the states to their trim values in about 3 
seconds. 

This example essentially demonstrates the effect of the second-order terms in the aircraft 
dynamics which are ignored in the linearized model. Just as in stability analysis, the influ- 
ence of higher-order terms becomes crucial when the first-order effects diminish. Since the 
symmetric thrust modulation of the PC A generates pitch moment variation and differential 
thrust change produces yaw moment variation, but no direct roll control is provided, the 
PC A would in this aspect resemble a crippled aircraft that loses the use of its aileron which 
provides the primary roll control. In this sense, the above exercise illustrates that when the 
cross-coupling nonlinearity in the system is appropriately taken into account in the controller 
design, it may actually renders the PC A more controllable. 

5. PCA Tuck-Mode Controller for the C-17 Aircraft 
5.1 Tuck-Mode of the C-17 Aircraft 

The C-17 is a heavy military transport aircraft and has four engines. A study was 
conducted in 1996 at the McDonnell Douglas Aircraft Company and NASA Dry den on 
engine-only flight control for the C-17. The study found that in most flight conditions the 
PCA controller can control the aircraft reasonably well. But at one particular condition 
specified by Mach 0.7, h = 30,000 ft, weight=576,000 lb, and zero flap, the longitudinal 
PCA controller failed to control the climb rate (or equivalently the flight path angle) to 
track h-commands. A closer examination of the linearized PCA model reveals that the C-17 
longitudinal dynamics at that condition have the so-called “tuck mode” [17]. This is the 
case when the phugoid mode degenerates into two real first-order modes, one of which is 
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unstable. Furthermore, the dynamics from the engine thrust to flight path angle at that 
condition are nonminimum-phase. In addition, it turns out that a singularity exists in this 
case in the system dynamics so that the conventional nonzero set-point tracking control 
design is not applicable (see next section). All these, in conjunction with the usual limited 
control effectiveness of thrust on aircraft attitude, render the controller design problem 
a challenging one. Part of the effort in this research was devoted to this problem. A 
satisfactory PC A controller was designed to provide accurate longitudinal control. Perhaps 
more importantly, the methodology used in this exercise, which combines the well- developed 
state-space techniques and state-of-the-art control design tools, may prove to be equally 
valuable in other applications. 

5.2 Optimal PID Control Laws 

In this section we briefly review an approach for design of optimal PID, full-state feedback 
controllers by an LQR method [18]. Given a linear system 


Ax + Bu 

(25) 

Cx 

(26) 


where x £ R n , u € /T, y € R m , and m < r < n. A nonzero set point z G R m is given for 
the system. This means that there exist a pair of constant steady-state {a: s ,u s } such that 
and u s satisfy Ax s + Bu s = 0 and y s = Cx s — z. To design a PID control law for y 
to asymptotically track z t Parker proposed the following technique [18]: Define a new state 
vector £ = ( y T r) T ) T = Tx, where y € R n ~ m and y = Lx for some (n — m) x n matrix L. 
The selection of L should make the transformation matrix 


T — 


^ c ^ 
\ L ) 


(27) 


nonsingular (which automatically requires C to be of full row rank). Then, 


i = TAT- 1 ^ + TBu = F£ + Gu 


(28) 
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Next, introduce the augmented state 9 = (9( 9\ @3 ) T £ R n+m , where 

0i = e = z - y, 0 2 = e, 0 3 = rj (29) 

and the new control v = u. The augmented state equations are then 

0 = $0 + Tv (30) 

where the matrices $ and T can be easily obtained from the definition of 9 and (28). Now 
the following LQR problem is solved where the performance index 

1 r°° rr, 

J — - / ( 9 t Q9 + v T Rv)dt (31) 

2 Jo 

is minimized. The optimal solution is given by 

u = v = —R~ x T t P0 = + /T 2 <9 2 + /T 3 ^3 (32) 

where P is the solution of the standard algebraic Riccati equation associated with (30) and 
(31). Integrating (32) once yields a control law 

u(t) — I< 1 [ e(r )dr + K 2 e(t) + I< 3 r)(t) (33) 

If 77 is selected in such a way that it contains the first-order derivative of y, this is a PID 
feedback control law. 

Remarks : 

1 . Although not mentioned explicitly in Ref. [18], the controllability of the augmented 
system (30) is obviously required for the unique positive definite solution of the asso- 
ciated algebraic Riccati equation. 

2 . Again, although not discussed in Ref. [18], the closed-loop system stability under 
this control law can be intuitively seen from how the control law is derived. The 
minimization of the performance index (31) ensures that 9 — V 0. By the definition of 
9 , this leads to e — * 0 (y = z) and y — > rj s where y s is some constant vector. Since 
rj — Lx and y = Cx, the steady-state value of x is obtained from x s — T~ l (z T yJ) T , 
which is necessarily finite. 
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5.3 Longitudinal Controller Design Procedure 

The complete 10-state, 6-DOF system equation for the C-17 includes some coupling be- 
tween the longitudinal and lateral/directional dynamics. To simplify the preliminary design, 
we shall at first ignore these coupling effects and the engine dynamics (which now are the 
actuator dynamics). Therefore the first-cut design of the controller is based on the 5-state 
longitudinal dynamics 

•Elon — ^-lon^lon + BlonU (34) 

y = Ci on xi on (35) 

with xi on = (0 q a v h) T and y = h = V 0 (0 — a), where the standard notation is used, i.e., 9 

denotes the pitch angle, q the body pitch rate, a the angle of attack, v the velocity, and h 

the altitude. The control is the engine throttle input. 

Aside from the tuck-mode and the nonminimum-phase behavior, this model also poses 
a special challenge for the tracking controller design: the conventional optimal nonzero set 
point tracking control law [19] does not exist because the matrix Ci on AY^ n Bi on ^ the inverse 
of which is required for the control law, is singular. In addition, it is desired to include an 
integrator in the controller to enhance the tracking performance, particularly in the pres- 
ence of possible system parameter uncertainties. Finally, the controller should use available 
measured feedback signals instead of requiring state feedback. 

The controller design was accomplished in two steps. First, the approach introduced 
in the preceding section was applied to the simplified longitudinal dynamics (34). The 
problem is posed as a tracking problem in which a constant h com is to be tracked by the 
output y = h. It should be pointed out that because of the dependence of the longitudinal 
dynamics on the altitude h, h = constant is not exactly a set (trim) point of the C-17 
with engine-only control. Also, although the system (34) is not strictly uncontrollable, its 
controllability matrix is nearly singular. This renders the solution of the algebraic Riccati 
equation inaccurate. When combined with the ignored coupling between the longitudinal 
and lateral dynamics and the engine dynamics, the resulting PID controller will not be 
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able to exactly control h to track h com . Nonetheless, the method described in the preceding 
section is systematic and easy to use. The outcome then serves as a starting point for further 
iterations in the controller design. 

Next, the controller obtained in the first step was applied to the complete 10-state sys- 
tem model with all the first-order engine dynamics and control rate-limiters included. The 
controller parameters were then further tuned by using the MATLAB Nonlinear Control 
Design (NCD) toolbox. The NCD toolbox allows the user to specify the desired shape of 
the response of a system with a graphic interface. The user-specified parameters of the 
system (which can be nonlinear) are then repeatedly adjusted by a sequential quadratic pro- 
gramming algorithm until the actual response of the system matches the desired one. Step 
response was used in the C-17 PC A controller design. As in any method based on numerical 
optimization, the initial values of the optimization parameters can be very important to the 
final success of the method. This is where the baseline design obtained in the first step fits 
in. 

Finally, the controller is switched from a state feedback controller obtained above to an 
output-feedback controller. At this point, the feedback signals to the controller, are xi on and 
the controlled variable h , and the control law takes the form of 

ll — k p (hcom T [ (hcom A)g?T R- x X( on (36) 

JO 

where K x £ R 5 . Among them, the measurement of the angle-of-attack a is usually quite 
noisy. So the normal acceleration a z (ft/s 2 ) measured at the Inertial Reference Unit (IRU) 
was used in place of a. The reason for using a z is that a z is approximately h = d(h — h com )/dt. 
Hence it provides the derivative information of the controlled variable h in the control law 
(the D-term in the PID control law). Let w = (6 q v h a z ) T . Then a linear transformation 
w = Hxi on exists, where the last row of the H matrix comes from the linearized relationship 
between a z and Xi on . The new gain vector with w as the feedback signal is simply K w = 
H~ l K x . Finally, the dependence of the control law on h stems from the weak dependence 
of the longitudinal dynamics on h. For the purpose of controlling h the A-feedback is not 
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necessary or desired. For instance, for an extended climb/descent, the magnitude of the h- 
feedback term will become larger and larger, and could eventually saturate the controller. To 
eliminate the h-feedback requirement, the corresponding component in K w was set to zero. 
The remaining gains, including k p and &*, were then fine tuned by using the NCD toolbox 
again to compensate for this final adjustment. The final output-feedback PC A controller 
takes the form of 

u — k p (Ja corn h) k{ f (hcom h)dr -f- kgO T kqQ T k v v T k az a z (37) 

J o 

Figure 9 depicts the controller configuration. Figure 10 show the response of the climb 
rate h to a 5 ft/sec command input. Both the transient and steady-state responses appear to 
be quite good. Notice the initial dip in the h history, which indicates a nonminimum-phase 
system. Figure 11 plots the corresponding engine command, which is not excessive. 
Remarks: 

1. The principle of the above design procedure can be applied to other cases. Standard 
LQR-based control design methods are systematic, easy to apply, and can yield very 
satisfactory performance. But they inevitably produce full-state feedback control laws 
which may not be realistic. On the other hand, the design of output-feedback con- 
trollers heavily depends on trial and error, and laborious numerical iterative processes 
are usually involved [20]. Empirical methods also exist for tuning of classical PID 
controllers. But for MIMO systems their applicability and effectiveness are severely 
limited. The design approach used here suggests that a different route may be taken. In 
this approach an appropriate LQR-based method is first applied, possibly to a simpli- 
fied but essential part of the complete dynamics. Once a baseline controller is obtained, 
it is added to the complete system model with actuator dynamics included. Then a 
state-of-the-art interactive graphic tools for control design, such as the NCD toolbox, 
is employed to tune the controller parameters. At this stage, appropriate measured 
outputs replace the unavailable state variables, and nonessential feedback paths are 
removed. The finally obtained controller should meet the design objectives. 
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2. It is recognized that this procedure does not automatically guarantee a success, and 
engineering judgments are still important in the process. But past experience [21] 
and the current exercise have indicated that this approach can significantly reduce the 
effort and time consumed in the intensive trial- and-error design iterations, and the 
controllers so obtained are usually quite satisfactory. 

6. Conclusions 

In this research the nonlinear phenomena observed in previous PC A flight testing were 
examined and analyzed. The predominant nonlinearities come from the engine dynamics, 
interactions between propulsion system and the airframe, and the dynamic cross-coupling 
effects between longitudinal and lateral-directional motions. It is concluded that if in the 
design of the PC A controller the essential aspects of these nonlinearities are incorporated, it 
is likely that the performance of the PCA controller can be markably enhanced to give more 
accurate flight control of the aircraft. A nonlinear predictive control method is introduced 
as a potential approach for design of such a nonlinear PCA controller. The problem of 
controlling a crippled aircraft, which from a control point of view bears some similarity with 
PCA, is solved to demonstrate that higher level of control effectiveness can be achieved when 
the system nonlinearities are included in the controller. 

As a part of the effort in this research and in conjunction with a PCA project conducted 
in the summer of 1996 at Dryden, a PCA longitudinal controller design for the C-17 trans- 
port aircraft was designed. In this process, although the PCA model was still linear, several 
attempts by using conventional linear control methods failed because of the particular dy- 
namic behavior of the aircraft at that flight condition. A methodology was developed that 
combines state-space based optimal PID control law design appraoch and interactive numer- 
ical design approach facilitated by the Nonlinear Control Toolbox in MATLAB. A successful 
controller was obtained for the C-17, and the methodology itself may prove to be a valuable 
addition to the existing design techniques. 
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Fig. 4: State histories of the aircraft for the longitudinal maneuver 



Fig. 5: Control surface defections for the longitudinal maneuver 
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Fig. 6: State histories when the linear controller applied to the nonlinear aircraft dynamics 
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Fig. 8: Control surface deflections for the turing maneuver 
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Fig. 9: C-17 Longitudinal PCA Controller 
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Fig. 10: Step response of the climb rate of the C-17 (command 
input = 5 ft/sec) 
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Fig. 11: C-17 engine throttle input corresponding to a commanded 
climb rate input of 5 ft /sec 



